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Abstract 

We consider a symmetric exclusion process on a discrete interval of S 
points with various boundary conditions at the endpoints. We study the 
asymptotic decay of correlations as S ^ oo. The main result is asymptotic 
independence of a stationary distribution whem the points are far away 
from each other. We develop a new recurrent probabilistic approach which 
is an alternative to Derrida's algebraic technique. 

1 Introduction 

Exclusion processes are a long-standing popular subject in both mathematics 
and theoretical physics. They are a particular case of processes with local inter- 
action on a lattice. The interest to them is mainly due to the fact that exclusion 
processes are the simplest nontrivial model for collisions in a multiparticle sys- 
tem. Using them, heat conduction [T], viscosity |2l|3j, quantum ferromagnet [4], 
nonequilibrium process [5], etc., models are constructed. Even irrespective of 
physics, these processes are natural probabilistic objects. The first systematic 
exposition of the relevant theory was given in the well-known monograph [S]. 

The simplest exclusion process is a symmetric exclusion process on the in- 
teger lattice Z. This process is well known to have a continuum of (spatially) 
uniform invariant measures, which are mixtures of Bernoulli measures. The 
same holds for an exclusion process on a finite segment of the lattice if bound- 
ary conditions are empty. For other boundary conditions, as a rule, it is easy 
to show that there are no invariant Bernoulli measures. Using a matrix method 
(ansatz; see [5l[7]), one can obtain an explicit form of correlation functions for 
such a process. This would imply that the invariant measure is asymptotically 
Bernoullian. It should be noted that this powerful algebraic method, similar 
to the famous Bethe ansatz, is rather cumbersome, is not always mathemati- 
cally well-grounded, and has substantial restrictions on its range of applicabil- 
ity. Namely, we are unaware of any application where jumps of particles can be 
longer than 1. It is also important that the probabilistic nature of the method 
is absolutely unclear. 
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In the present paper, we propose another — simple and natural from the prob- 
ability theory viewpoint — approach, which can be extended to jumps of lengths 
greater than 1, as well as to other boundary conditions. Being quite differ- 
ent from the Bethc-Derrida methods, our approach has something in common 
with them, namely, a certain recursive procedure. Here we demonstrate the 
idea of our approach in the simplest situation. Various generalizations will be 
considered in further publications. 

2 Problem Statement and the Result 

We consider a simple symmetric exclusion process on an interval Is ~ {0, 1, . . . , S', S+ 
1} of a ID lattice. The state space for this continuous-time finite Markov chain 



is the set {0, 1}^^. Jumps arc defined as follows. For any time interval [t, t + dt] 
and any pair s,s + l, where s — 0,1, . . . , S, the nodes s, s -|- 1 interchange their 
values with probability Xdt (independently for different s), so that 



This defines a generator of a Markov process that we call a process with empty 
boundary conditions. However, we shall consider boundary conditions: 



More precisely, this means that ^o(O) = and interchange between the points 
0, 1 G Is reduces to the following: the pair ^f(O) = 0, Ct(l) = 1 becomes 
Ct+aft(0) = 0, £,t+dt{^) = with probability Xdt. The same for the points 
S,S + 1 e Is: ^o{S + 1) = 1 and the pair Ci{S) = 0, Ct(5' + 1) = 1 becomes 
^t+dt(S) = 1, ^t+dtiS + 1) = 1 with probability Xdt. 

Denote by 7r(no, rii, . . . , ns+i) = ■7r^^\no, ni, . . . , ns+i) the stationary mea- 
sure of this Markov chain. We are interested in the functions 

m^'^^(a;i,a;2,...,a;fe) = n^^^n^^ = 1, n^,^ = 1, . . . , n^^ = l), xi < ... < Xk- 

For the empty boundary conditions, the process is reducible (because of the 
conservation of the number of particles, i.e., the number of ones), and any 
stationary measure is a mixture of Bernoulli measures, i.e., measures that are 
products of S+2 equal measures. As is easy to show, in the case of our boundary 
conditions the process is irreducible (the number of particles is not preserved 
now), and its invariant measure is not a product of independent measures. Our 
goal is to prove the following correlation decay property for this process. 

Theorem 1 For allO < x < S + 1, we have 



(6(0), 6(1), •••,6(5 + 1)), 6(^) = o,i 



6+dt(s) = 6(s + i), 



6+dt(.s + l) = 6(s). 



6(0) = 0, 



6(<5 + i) = i. 



nil (^) 



.r 



5+1' 
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For any < a;i = 
as S ^ oo, we have 



xi{S) < X2 — X2{S) < S + 1 such that 



S+i 



ak, fc = 1, 2 



3 Proof 

The proof consists of two parts. The first part, reduction to a dual process 
(which is easier to study), is rather standard but cumbersome. The second 
part, analysis of the dual process, is original. 

3.1 Moment Closeness 

First we obtain equations for moments of the process fj: 

k 

mfc(i; {xi, . . . = E]^ff(a;i), xi, . . . ,Xk & Is, fc = l,2,.... 

1=1 

For a process with values and 1, these moments uniquely determine marginal 
distributions 



and, of course, vice versa. These equations possess the property of moment or 
marginal closeness [8], which makes it possible to construct a convenient dual 
process. In what follows, we assume for simplicity that A = 1. 

We introduce the notion of a cluster. Let a set Q C N be given. A subset 
(fci < ... < km) — K Q is called a cluster if fc^+i — fc^ = 1 for all i = 1, . . . , m— 1 
and, moreover, {fci — 1} ^ Q and {km + 1} ^ Q. In other words, a cluster is a 
maximal chain of successive elements of Q. 

Lemma 2 Consider a set X — (I'l, . . . , x^) such that < xi < . . . < Xk < 

p 

15*+ 1. Let X = Kj, i.e. X is a union of clusters Kj — (^x\''^\ . . . ,x'f \ ). 
i=i 

Then for any X we have the equation 

— mfc(i, {xi, . . .,Xk}) = ^mfe(i, {xi, . . .,x\l^ - 1, ■ ■ ■,Xk}) 
p 

+ ^'rnk{t, {xi, . . • + 1, . . . - 2pmk{t, {xi, . . . ,Xfc}) 

with boundary conditions 

mfc(t, {0, X2...,Xfe}) =0, 
mk{t, {xi . . . ,Xfc_i,5 + 1}) = mfc_i(i, {xi . . . ,Xfc_i}). 
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Proof. First wc derive the equation for a one-point function. The moment 
closeness property is based on the Unearity of the conditional probability 

Pr{^t+dt{x) = 1 I ^(x-l) = a,^t{x) = /3,Ux+l) = -f) = adt+P{l-2dt)+jdt. 

By the total probability formula, 

1 

Pr{it+dt{x) = l)= Yl {dt{a + ^-2/3)+/3)Pr{^t{x-l) = a,Ux) = P,Ux + l)=-f) 

a,/3,7=0 

= dt(^Pr{^t{x - 1) = 1) + Pr{^t{x + 1) = 1) - 2Pr{^t{x) = l)) + Pr{Ux) = l), 

whence we get 

^mi(t; {x}) = mi{t; {x + 1}) + mi(t; {x - 1}) - 2mi(t; {x}). 

Now we pass to correlation functions of an arbitrary order k. Let X consist 
of clusters of length 1 only. Since a change during time dt can only occur at one 
of the points Xi, similar computations for each of them yield 

d ^ 

■j^mk{t; {xi,...,Xk}) = ^(rnkit; {xi,...,Xi + l,...,Xk}) 

i—l 

+ mk{t] {xi, . . . ,Xi + 1, . . . ,Xfc})^ - 2kmk{t; {xi, . . . 

Now consider the situation with clusters of larger lengths. As above, it suffices 
to consider the case of only one such cluster. Thus, let Kj = (x^^^ , . . . , a; -^^^^.^ ) be 
the jth cluster in the partition of X. Changes can only occur at the endpoints 
of the cluster. Let us write the general form of the conditional probability in 
this case, omitting the index j: 

Pr(^^t+dt{xi^) = 1, • • • , ^t+dt{xi^) = 1 I ^t{xi^ - 1) = a, 
^t{xn) = /3, = 7, + 1) = ^) = ajdt + 0j{l -2dt) + pddt. 

Of course, the conditional probability is taken provided that ^tixi/) = 1, Vj = 
2, . . . ,q — 1. Then we obtain 

Pr{£,t+dtixtJ = 1, . . . , S,t+dtixiJ = 1) 
1 

= Yl {<^ldt + (3^{l-2dt)+(35dt)Pr(^^t{xn-l) = a,^t{xn) = P, 

a,/3,7,(5=0 

^t{xi, + 1) = 1, . . . , - 1) = 1, 6K) = 7, + 1) = '^j 

= dt(mg{t; {xi^ -l,...,XiJ)+mq{t; {xi^,...,Xi^ +1}) 

- 2mg{t; {x^^ ,...,XiJ)] + mq{t; {xi^ ,...,x^J), 
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which implies that for an individual cluster we have the equation 

^mq{t■,{x,,,...,X^J) =mq{t;{xi^ - 1,...,X^J) +mq{t■,{x^,,...,X^^ + 1}) 

-2m,(i;{x,i,...,x,J). 

Taking the sum over clusters, we obtain the lemma. For convenience, we write 
the right-hand side of the equation as follows: 

d 

— mfc(t; {xi, . . .,Xk}) ^^Al^rukit; {xi, . . .,Xk})- 
3.2 Dual Process 

Denote by 2^^ the set of (finite) subsets of Is, including the empty set 0. Let 
i^t, t > 0) be the exclusion process defined above. A process {At, t > 0) with 
values in 2^^ is said to be dual to the process (^t, t > 0) if for alH > we have 

E n ^*(^) = ^ n ^oix). 

xGAq x^At 

In [6l Section VIII. 1, Theorem 1.1] it is proved that a symmetric exclusion pro- 
cess is self-dual. Thus, for instance, analysis of two-point correlation functions 
can be reduced to consideration of two particles performing an exclusion walk. 
To extend this result to our boundary conditions, consider the above equations. 
They lead to the following statement. 

Lemma 3 A process dual to (^t, t > 0) is an exclusion process [At, t > 0) for 
a finite number of particles with the following modification: 

INI 

1. If one of the particles of At touches the boundary 0, the whole configuration 
dies (we also say that all particles reach the absorbing state 0): 

a := sup{s > : n {0} 7^ 0} =^ A„ = 0; 

2. Any particle that reaches the boundary + 1 sticks to it. After that, the 
other particles walk independently of the stuck particle, whose coordinate 
remains to be S + 1. 

Note that after one of the particles reaches the right-hand boundary, other 
particles either hit or reach the point 5*+ 1. Therefore, with probability 1, the 
Markov process thus defined comes in the course of time to one of the absorbing 
states: or {5 -f- 1, . . . , S* + 1}. 

As a consequence, we get the following statement. 

Proposition 4 We have 

mf\A)^ ^ Pr(£), \A\^k, 
eeL{A) 
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where L{A) is the set of all trajectories of the dual process that start in A and 
do not hit the boundary 0. 

Proof. Let the process rit{x) = (rjl {xi), . . . ,r]^{xk)) S be an exclusive 
random walk of A; particles in the segment [0, 5+1], where x = {xi, . . . ,Xk) G 
is the initial disposition of particles. To prove the duality, it suffices to show 
k 

that the function E Yl ^oiVti^i)) satisfies the same differential equations as the 

»=1 k 

correlation function E fl ^t{xi)- By the total probability formula, we have 

i=l 

fe 

Ell^oiviixi)) = Yl P^i^oik) = l,...,Co(ifc) = l)Pr{'nl{x,) = i,,...,r,'^{xk)=ik). 

i=l ii<...<ik 

Denote 

Pt{ii, ...,ik)-= Pr{ril{xi) = ii, . . . , Tni^{xk) = in)- 
This function satisfies the same equation as the moments, i.e., 

d 

—Pt{ii,...,in) = ^A?.Pt(ii,...,ife). 
Indeed, for a single particle we have 

i+l 

Pt+dt{i) = X] Pr{r]t+dt{x) =i\r]t{x)= j)Pr{r]t{x) = j) 

= dt(Pr{rit{x) =i-l)+ Pr{r]tix) = i + l)) + (1 - 2 dt)^{r]t{x) = i). 

Hence, 

^PtW =Ptii-l) + Pt{i + 1) - 2Pt(i). 
Similar arguments also apply to the case of several particles. Thus, 

k 

j^Efl^Mi^i)) = J2 Pr{U^l) = h .■.,a^k) = l)j^Pt{ti,-.-,tk) 

i=l ii<...<ik 

ii<...<ife lj=l ) 

= E |E^I^H?o(zi) = 1, ?o(»fc) = l)|pt(»i,...,*fe) 

ii<...<Zfe li=l ) 

k k 

j=i i=i 
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Clearly, boundary conditions are the same as for the original exclusion process: 



fe-1 



i=l 

fe-1 fe-1 

E n ^o{vl{xi))^o{vl{S + 1)) = n ^oivlixi)). 

i=l i=l 

Therefore, the process that we have constructed is indeed dual to the original 

exclusion process. Summing up and taking into account the fact that all particles 
finally get into one of the two absorbing states, we obtain 

k k 

mk{oo;{xi,...,Xk}) = -B J|Coo(a;j) ^ EY[^o{vlo{xt)) 

i=l i=l 

= Pr{v'M = S+l, vL{xk)=S+l)=P^\ 

where is the probability that all particles in the dual process reach the 
absorbing state 5 + 1. 

3.3 Analysis of the Dual Process 

Consider the exclusion process r]t = r]l°°^ = {xt,yt) for two particles on the 
segment [0, 5+1]. At the initial time instant, the particles are at points < 
Xo{S) < ya{S) < 5+1, respectively, such that — >s^oo « and -^s^oa 
j3. It is also convenient to assume that the distance between them is at least 2. 
Clearly, this does not lose generality. Let us show that 

Pg^ ^s^oo a/3, 

whence the theorem follows. 

For an exclusion process ryt, we introduce time intervals {TkjTj,) when par- 
ticles are at distance 1 from each other, assuming that the rest of the time 
the distances are greater than 1. We call the fcth meeting moment. Using 
these time moments, for k = 0,1,... we define processes ■q^^\t) (which are 
non-Markovian if fc > 0!). In each of the processes, two particles with the same 
initial disposition are walking. 

In the process ri^^\t), particles walk independently (simple symmetric ran- 
dom walk) and do not see each other. In the process ri^''^t), the exclusion 
process works up to time moment Tj.. After that the particles become at dis- 
tance 2 from each other, and from this point they start walking independently 
and do not see each other. 

Let Pfe be the probability that in the process rj^''^ {t) both particles come to 
the point 5+1 before one of them hits the origin. Since the number of particle 
meetings before absorption is finite with probability 1, it sufiices to prove that 
for large 5 the probabilities Pk differ little from Pq. 
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First of all, it is easy to prove that asymptotically as 5 ^ oo we have 



Our main method is recurrence relation between the probabilities Pk- 

First we compare Pi and Pq. Define a random variable t(x, y), the moment 
of the first meeting of the particles provided that they started from the points 
X and y. It is clear that for all processes rj^'^^t), k > 0, the distribution of this 
variable is the same. Let Xr and = a;,- + l be points where the first and second 
particles, respectively, are at the first meeting moment, and let ^oi^r = n) be 
the probability that the first particle is at point n at this meeting moment. Let 
Pk{x, y) be the probability for two particles in the A;th problem to reach 5* + 1 

if at the initial time moment they are at points x and y. Note also that the 

s 

event IJ (.x^ = n) comprises also the outcome that the two particles reach 5+1 

?i=i 

earlier than die, since in this case the distance between them would necessarily 
be equal to one, at least at the points {S, S + 1). 
Then, by the total probability formula, 

1 S 
A = 2 I] Poi^r = n)(Po(n,n + 2) + Po(n - 1, n + 1)) + ^-^Po{xr = S). 

n=l 

We transform the sum in parentheses, first rewriting each term: 

n(n + 2) n(n +1) n 



Po(n, n + 2) - Po(n, n + 1) = 
Po(n-l,n+l)-Po(n,n+l) 



(5 + 1)2 (5 + 1)2 (5+1)2' 

(n-l)(n+l) n(n + l) n- 



(5+1)2 (5+1)2 (5+1)2' 
and therefore 

Po(n, n + 2) + Po(n - 1, n + 1) = 2Po(n, n + 1) - ^ . 

(,o + ij 

Thus, we finally obtain an expression for Pi via Pq: 

Pi = Po{xr = n) (^Po(n, n + 1) - + -^P^i^r = S) 



= Pf 



S , S-l 



^ Poixr = n)Po(n,n+ 1) - -^Tg^^-^^ Pq{xt = n) 

n=l ^ ' n=l 



2(5+1)2- 



To getgtlje last equality, we again use the total probability formula. Note that 
Ci = Yli Po{xt = n) is a positive constant less than one, since this sum is the 



n=l 
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probability that the particles meet at least once before one of them is absorbed 
at or 5* + 1. 

Now, taking the obtained result into account, we make similar reasoning for 
the probability P2: 

1 S 
P2 = 2 E ^o(^- = n){Pi{n,n + 2)+P^{n-l,n+ 1)) + j:^Po{xr = S) 

1 ■^"^ / 

= ;r V Po{xr = n) Po{n, n + 2) + Po{n - 1, n + 1) 

^ n=l V 

- 2{S + ir (^1 " + 2) + Ci (n - 1, n + 1))) + {xr = S) , 

where 

s-i 

Ci(j,j) = 51 ^o(a;r(i,j) = n). 

n=l 

Denote 



1 ^"^ 

= 2 ^o(a;r = n) (Ci (n, n + 2) + Ci(n - 1, n + 1)) ; 



71=1 

then 



P2=Pl 



2(5+1)2 

Taking into account that C\ < 1, we roughly estimate C2: 

^ 5-1 s-i 
C2 = -EPo(x^ =n)(Ci(n,n + 2) + Ci(n-l,n+l)) < ^%(x^ = n) < 1. 

71=1 71=1 

Analogously, a similar relation can be written for any fc: 

* 2(5+1)2' 

where 

1 

Cfe(«)i) = 2 X] ='^)(Cfe-i(">'^ + 2) + Cfe_i(n-l,n+l)), 

71=1 

Cfc :=Cfc(a,/3), Cfc < 1. 

We rewrite this as follows: 

k 



^'=-^°~2(5TTF§^^' 



Passing to the limit as A; — )• oo, we obtain 



^ ' i=l 

Thus, we have an expression for the difference of probabilities 

^ oo 
^ ' i=l 

It remains to estimate this difference and show that it tends to zero as S ^ oo. 

By the construction, Ck is the probability that two independently walking 
particles meet at least k times before one of them is absorbed. To estimate this 
quantity, we introduce an auxiliary model. 

Consider a process Xt € [0, S], symmetric random walk of a particle on the 
segment [0, S] with Xq = 1. Let 7^ be the probability for a particle to come to 
zero k times before reaching S. We claim that is not greater than 7^,. Indeed, 
Xt + 1 can be interpreted as the distance between two particles in our original 
problem at time instant t. The condition Xq = 1 means that at the initial time 
moment we place the particles at the smallest distance before their meeting, 
i.e., at distance 2. If Xt = S, one of the particles is absorbed for sure, since 
the distance is equal to 5+ 1. Therefore, if two independently walking particles 
meet k times before one of them is absorbed, then the auxiliary particle Xt hits 
zero k times before it comes to S. Thus, the probability to hit zero k times in 
the auxiliary model is (at least) not less than the probability to meet k times 
in the original problem. 

Let us calculate 7^. It is clear that 71 = If a particle comes to 0, it 

makes the next step to the right, to the point 1, and everything starts over 
again; hence, 

5-1 

7fe = 7fc-i ■ 

Finally, we obtain 

5-1 

Ik 



5 

Now we can again estimate the difference between Pq and Ptx 
P0-Poo= ^ E ^ 2(5 + 1)2 I^^' 

1 ^fs-iY 5-1 



2('S + l)^^V S J 2(5+1)2 2(5 + 1) (5+1)2' 

Therefore, [b]|Poo — -Pol = O^^q-j-^ as 5 ^ 00; since \b]Pk -^k^oo Poo for any 
5, this completes the proof. 
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